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ABSTRACT 

We present high signal-to-noise optical spectra of two luminous super star clusters in 
the star burst galaxy M82. The data for cluster F and the nearby, highly reddened 
cluster L were obtained with the William Hcrschel Telescope (WHT) at a resolution 
of 1.6 A. The blue spectrum (3250-5540 A) of cluster F shows features typical of mid- 
B stars. The red spectra (5730-8790 A) of clusters F and L show the Call triplet 
and numerous F and G-type absorption features. Strong Call and Nal interstellar 
absorption lines arising in M82 are also detected, and the A6283 diffuse interstellar 
band appears to be present. The quality of the WHT spectra allows us to considerably 
improve previous age estimates for cluster F. By comparing the blue spectrum with 
theoretical model cluster spectra using the PEGASE spectral synthesis code (Fioc 
& Rocca-Volmerange 1997), we derive an age of 60 ± 20 Myr. The strength of the 
Call triplet is also in accord with this age. Cluster L appears to have a similar age, 
although this is much less certain. The measured radial velocities for the two clusters 
differ substantially, indicating that they are located in different regions of the M82 
disk. Cluster F appears to be deep in M82, slightly beyond the main star burst region 
while the highly obscured cluster L lies near the outer edges of the disk. We derive 
an absolute magnitude My = —16.5 for F indicating that it is an extremely massive 
cluster. The presence of such a luminous super star cluster suggests that the M82 
starburst experienced an episode of intense star formation ss 60 Myr ago. 

Key words: galaxies: evolution - galaxies: individual (M82) - galaxies: starburst - 
galaxies: star clusters - galaxies: stellar content 



1 INTRODUCTION 

M82 is one of the nearest examples of a spectacular starburst 
galaxy. At a distance of 3.6 Mpc (see Freedman et al. 1994) 
its far infrared luminosity is about 3.2x10 (Colbert et 
al. 1998), and the many effects of the starburst, including 
extreme optical surface brightness and giant outflows, are 
immediately apparent (e.g., Lynds & Sandage 1963). M82 
has naturally played a special role in studies of intense star 
formation on large scales (e.g., Telesco 1988). Among the re- 
sults from ground-based optical and infrared investigations 
of M82 was the recognition that dense, luminous star clus- 
ters are common in M82 (O'Connell & Mangano 1978; here- 
after OM78). High angular resolution optical imaging with 
the Wide Field Planetary Camera on the Hubble Space Tele- 
scope (HST) revealed more than 100 clusters located in the 
optically visible skin of the M82 starburst region (O'Connell 
et al. 1995). Infrared imaging penetrates the dense gas and 
dust in the starburst and finds populations of compact, lu- 
minous star clusters in the actively star-forming regions of 
the galaxy (Satyapal et al. 1995, 1997). Star clusters are 



evidently a featured product of the intense star formation 
now occurring within M82. This is consistent with the more 
extensive results from the HST showing compact, luminous 
'super star clusters' (SSCs) to be frequently associated with 
starbursts under a wide range of conditions (e.g., Whitmore 
& Schweizer 1993; Whitmore et al. 1995; Conti & Vacca 
1994; O'Connell, Gallagher & Hunter 1994; Meurer et al. 
1995; Schweizer et al. 1996; see Ho 1997 for a review). 

The presence of super star clusters in starburst galaxies 
provides us with a useful tool for exploring some aspects of at 
least their recent star formation histories, while also yielding 
some information about chemical abundances in stars. These 
measurements are feasible because star clusters closely ap- 
proximate a coeval, single metallicity simple stellar popula- 
tion (SSP). An SSP is the simplest class of stellar system 
to model, and one in which age and metallicity can be ex- 
tracted from the integrated spectrum (see Bruzual & Char- 
lot 1993, Bruzual 1996). In return for their increased com- 
plexity, SSCs offer several practical advantages as compared 
with observations of individual stars; most significantly they 
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can be 10-100 or more times as luminous as the brightest 
supergiant stars, allowing SSCs to be observed at greater 
distances or against more complex backgrounds. 

Achieving the potential of SSCs as astrophysical probes 
of starburst galaxies, however, has proven difficult. Star- 
bursts are often dusty, limiting the utility of optical and 
ultraviolet colors to yield ages of younger clusters (e.g., Wat- 
son et al. 1996). For intermediate age or older star clusters 
the well-known age-metallicity degeneracy becomes an issue 
(Searle, Wilkinson & Bagnuolo 1980). Optical and ultravi- 
olet spectra are available for a few blue SSCs, but initial 
studies were limited to qualitative discussions (e.g., Arp & 
Sandage 1985, Lamb et al. 1985). More recent observations 
have led to considerable progress in analyzing very young 
SSCs, where strengths and profiles of ultraviolet atomic res- 
onance lines are powerful diagnostics (e.g. Vacca et al. 1995, 
Gonzales Delgado et al. 1997). 

Optically luminous SSCs tend to be dominated in the 
blue by B- or A-type spectra (Schweizer & Seitzer 1993, Zepf 
et al. 1995; Brodie et al. 1998). As a result, ages are not 
readily determined from measurements of individual spec- 
tra features, especially when the Balmer jump is not avail- 
able (Diaz 1988). In addition, these studies have found that 
the observed strengths of the Balmer lines exceed those pre- 
dicted by the Bruzual & Chariot (1993) models. Brodie et 
al. (1998) find that a truncated IMF is needed to reproduce 
their observations of the galaxy NGC 1275. For the M82 
starburst, Satyapal et al. (1997) use a wealth of infra-red 
data to infer a typical age of 10 7 yr for the clusters within 
a radius of 250 pc of the nucleus. They also find an age dis- 
persion of 6 x 10 6 yr which is correlated with the projected 
distance of the clusters from the nucleus, and therefore sug- 
gest that the starburst is propagating outwards from the 
nucleus. 

In this paper we present a new spectroscopic study of 
the luminous SSC F in M82 (OM78) and the nearby, highly 
reddened cluster L (Kronberg, Pritchet & van den Bergh 
1972). Cluster F has a size of about 9x5 pc from observa- 
tions with the original Planetary Camera on HST and an es- 
timated luminosity of My » -14.5 (O'Connell et al. 1995), 
and is therefore one of the most luminous younger star clus- 
ters known. It is located 25 arc sec (= 440 pc) south-west 
of the 2.2 /jm nucleus and is considered to be in a less ob- 
scured region of the starburst (O'Connell et al. 1995). The 
goal of our program is to develop improved approaches for 
quantitatively matching spectra of SSCs to model predic- 
tions as a means to determine cluster ages. Observations 
for this project were made at the 4.2 m William Herschel 
Telescope (WHT) using the double beam spectrograph ISIS. 
This allowed us to obtain wide wavelength coverage with 
good spectral resolution, and thereby considerably improve 
the estimated age for cluster F. 

The observational procedures are reviewed in the next 
section. In Sect. 3 we describe the spectra, while Sect. 4 com- 
pares observed spectra with predictions from a population 
synthesis model. Results are discussed in Sect. 5, including 
implications for the duration of the M82 starburst and the 
quality of current models for SSPs with near solar metallic- 
ity. 



Table 1. Journal of WHT Observations. 

Blue Arm Red Arm 

Wavelength Exposure Wavelength Exposure 
Range Time Range Time 

(A) (s) (A) (s) 



3234- 


-4042 


2400 


5726- 


-6534 


1200 








6478- 


-7286 


1200 


3984- 


-4792 


1200 


7224- 


-8032 


1200 


4729- 


-5537 


2400 


7977- 


-8785 


2400 




30 



20 



10 -10 

Offset (arcsec) 



-20 



-30 



Figure 1. The raw spatial profiles for the red wavelength region 
(8100-8600 A; solid line) and the blue region (4100-4600 A; dotted 
line) normalised to the same exposure time. The two superclusters 
F and L are indicated. The x-axis represents the spatial offset 
with respect to the centre of SSC F and SSC L is 4" to the SW 
of F. The dashed line shows a representation of the polynomial 
fit used to subtract off the underlying galaxy background for the 
red spectrum. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of SSC F in M82 were obtained on 25/26 
March 1997 with the WHT on La Palma, Canary Islands. 
SSC F was observed at a position angle of 20° so that spectra 
of SSC L could also be obtained. The journal of observations 
is given in Table |l]. The double-beam spectrograph ISIS was 
used with Tektronix CCDs (1024 x 1024 24^m pixels) on 
both the blue and red arms. Two 600 line gratings were used 
to give a continuous wavelength coverage from 3234-5537 A 
and 5726-8785 A. A long slit of dimensions l".0 x 4'.0 was 
chosen to give a two-pixel spectral resolution of 1.6 A. The 
spatial scale was 0". 36 pixel -1 and an average seeing of l".2 
was measured over the period of the observations. 

The data were reduced using the PAMELA optimal ex- 
traction routines (Home 1986) within the FIGARO (Short- 
ridge et al. 1997) software package. Each frame was first 
bias-corrected, and those obtained at the reddest wavelength 
setting were divided by a normalised flat-field to remove the 
fringing. In Fig. ji], we show the spatial profile along the slit 
at representative red and blue wavelengths obtained by col- 
lapsing the raw images from 8100-8600 A and 4100-4600 A. 
It can be seen that SSC L is not detected in the blue re- 
gion, and that there is a substantial contribution from the 
galaxy background. This was subtracted by using polyno- 
mial fits to each pixel in the spatial direction, as shown 
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in Fig. |lj. The object data were then optimally extracted 
and wavelength calibrated using Cu-Ar comparison spectra 
taken before or after each integration. The red arm spectra 
were next corrected for telluric absorption features using a 
B star standard, and all spectra were extinction-corrected 
using the curve appropriate to La Palma. 

Finally the data were flux calibrated using the flux stan- 
dard Feige 34 (Oke 1990) and the individual wavelength 
ranges were merged. To accomplish this, small adjustments 
of 5-15% in flux level were required. The final spectra have 
a signal-to-noise ratio per pixel in the continuum of 30-40 
above the Balmer jump for SSC F and longward of 7000 A 
for SSC L. 



3 DESCRIPTION OF SPECTRA 
3.1 Continuum 

In Figs. 2 and 3 we show the flux calibrated spectra for SSC's 
F and L and enlargements over selected wavelength ranges 
to highlight various spectral features. We have performed 
synthetic photometry on the spectra to derive magnitudes 
and colours. For SSC F, we derive UBV magnitudes of 17.6, 
16.9 and 15.8 respectively giving (U — B) = 0.70 and (B — 
V) — 1.07. These colours are in excellent agreement with 
those of OM78 who derive V = 15.00, (U - B) = 0.69 and 
(B — V) — 1.04, using a 7". 5 aperture. From HST images, 
O'Connell et al. (1995) find V = 16.3 within a radius of 
0".5. For SSC L, we derive a monochromatic R magnitude 
of 16.7 compared to 15.6 for F. 



3.2 Absorption Lines 

Line identifications, heliocentric velocities and equivalent 
widths of the principle spectral features are given in Table 
2 for SSC's F and L. The equivalent widths were measured 
by defining continuum windows either side of the features 
and the errors given in Table 2 represent lcr uncertainties 
based on the signal-to-noise in the continuum windows. We 
also list equivalent widths based on the windows defined by 
Bica & Alloin (1986) which tend to be much wider than the 
actual spectral features present in our data. We now discuss 
each SSC in turn. 



3.2.1 M82 SSC F 

The blue spectrum is dominated by broad Balmer absorp- 
tion lines and a strong Balmer jump. Comparison of the 
equivalent widths of H9, H8 and H5 (least likely to be af- 
fected by nebular emission) with main sequence stars in the 
Jacoby, Hunter & Christian (1984) spectral library shows 
that the best match is for a mid-B spectral type. This clas- 
sification is confirmed by the detection of weak He I lines 
(e.g. A4922 to the red of H/3 in Fig. 2.) and Mgll A4481. 
We find that the strong and narrow Ca II K line arises from 
interstellar gas in M82 rather than stellar atmospheres since 
its velocity agrees with that of the Nal D lines and is very 
different to the stellar absorption features. 

In the red spectrum, the strongest absorption features 
present belong to the Ca II triplet and we also detect broad 



absorption from members of the Paschen series. The equiv- 
alent widths given in Table 2 for the Call triplet have been 
measured using the continuum and line windows of Diaz, 
Terlevich & Terlevich (1989). The combined strength of the 
two strongest members (Ca2+Ca3) is 8.5 A, uncorrected for 
Paschen line absorption. The mean velocity of the triplet is 
+35 kms -1 and probably represents the most accurate value 
for the radial velocity of SSC F since the Balmer lines are 
broad and contaminated by nebular emission. 

The overall spectrum redward of 5000 A has various fea- 
tures attributable to F and G stars (Andrillat, Jaschek & 
Jaschek 1995; Carquillat et al. 1997). In addition to the Call 
triplet, we detect weak absorption features due to CrI+Cul 
A5782, Fe I+Ba Il+Ca I A6497 as well as various Ti I and Fe I 
lines near the Call triplet. No molecular bands due to e.g. 
TiO are detected in our spectra. 



3.2.2 M82 SSC L 

The red spectrum (Fig. 3) clearly shows that the Ca II triplet 
is present and a weak Pa 12 absorption line is detected, 
suggesting that blue stars are present in this cluster. Weak 
low ionization absorption features arising from F and G stars 
are also detected and have similar strengths to those found in 
F. The strengths of the Call triplet lines are similar to those 
measured for SSC L with Ca2+Ca3 = 8.9 A. The radial 
velocities are, however, quite different with a mean for the 
Call triplet of + 131 kms -1 compared to +35 kms -1 for F. 
This large difference in radial velocity indicates that clusters 
L and F lie at different locations along our sight line through 
the disk of M82. 



3.3 Emission Lines 

For SSC F, we detect nebular emission lines of Ha, H/3, 
H 7 , [Oil] A3727,3729, [OIII] A4959, 5007, [OI] A6300, [Nil] 
A6548, 6584, and [S II] A6717, 6731. The H/3 and Ha + [Nil] 
emission line profiles are shown in Fig. 2 and heliocentric 
radial velocities and equivalent widths are given in Table 
2. For the H/3, Hy, [OI], [Oil], [OIII] lines, we derive a 
mean velocity of +181 ± 12 kms - which agrees well with 
the mean velocity of the Ca II K and Na I D interstellar lines 
of +189 ± Skins" 1 . The Ha, [Nil] and [S II] lines show a 
lower velocity of +114 ± 9 kms -1 in good agreement with 
that determined by OM78 of +121 ± 10 kms -1 for the same 
features. Inspection of the images shows that the spatial 
structure of these lines is more complicated than the [O II] 
or [O III] emission which is constant in intensity and veloc- 
ity across the region displayed in Fig 1. The emission for the 
strong Ha, [N II] and [S II] lines is clumpy with two veloc- 
ity components - the component seen in the other lines at 
+ 181 kms -1 and another component at +114kms _1 which 
lies between SSCs F and L. As a result of this structure, 
the strong +181 km s -1 component is over-subtracted from 
both SSC spectra and appears in absorption in Figs. 2 and 
3. 

For SSC L, we detect nebular emission lines of Ha, [N II] 
and [S II] and find compared to SSC F a lower mean helio- 
centric velocity of +69 ± 9 kms -1 . The equivalent widths of 
the nebular lines are similar for the two spectra. 
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Figure 2. The blue and red spectra of SSC F (smoothed with a = 1.0 A for presentation purposes). The lower three boxes (at the 
full spectral resolution) show enlargements covering the H/3, Ha+ [Nil], and Call triplet regions. The y-axis is in units of 10 — 15 
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Figure 3. The red spectrum of SSC L (smoothed with a = 1.0 A for presentation purposes). The lower three boxes (at the full spectral 
resolution) show enlargements covering the Nal interstellar absorption lines, Ha+ [Nil], and Call triplet regions. The y-axis is in units 
of 10 -15 ergs s _1 cm -2 A -1 



3.4 Diffuse Interstellar Bands 

Due to their moderate-to-high levels of interstellar dust ob- 
scuration, M82 clusters F and L provide excellent sight lines 
for measuring interstellar features in the spectrum of M82. 
We have already noted the presence of strong interstellar ab- 
sorption from the atomic species in the Call K and Nal D 
lines. Our spectra also have sufficient resolution and signal- 
to-noise ratios to allow sensitive searches for diffuse inter- 
stellar bands (DIBs). Major features of DIBs are reviewed 
by Herbig (1995). The strongest DIBs in the optical spec- 
tral region are the A5778-80, A6283, and A6533 features, and 
we therefore looked for these in our spectra. The A5778-80 
feature is coincident with the A5782 Cr I+Cu I blend which 
we positively identify in the spectra of F and L. The A6283 
DIB, is, however, probably present. We see a weak feature at 
6288 A (see Table 2) which is definitely real as it is occurs 



in the overlap region between spectra and is seen in both 
sets of data. Comparison with the spectral library of Jacoby 
et al. (1984) shows that it does not appear to be a stellar 
absorption feature. Likewise, comparison with the sky spec- 
trum and the atmospheric standard shows that it is not due 
to over-subtracted sky or a telluric feature. The strength of 
the A6288 absorption also correlates with reddening, being 
stronger in cluster L, as expected if this is a DIB. The pos- 
sibility therefore remains that the feature at A6288 could be 
a DIB arising from the interstellar medium of M82. 

The absorption strengths of DIBs are complicated. 
They are usually strongest in the boundaries of molecular 
clouds and tend to weaken in purely molecular regions as 
well as being absent from H II regions (e.g. Jenniskens et al. 
1994). This behavior is consistent with models where DIBs 
are produced by ionized large molecules (Sonnentrucker et 
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Table 2. Velocity and Equivalent Width (Wj) Measurements. 



Line A lab SSC F SSC L Comments 

Velocity Measured Window 1 Velocity Measured Window 1 

(A) (krns" 1 ) W A (A) W A (A) (krns" 1 ) W A (A) W A (A) 



[OH] 


3727 


+176 


-8.81 + 0.21 










Nebular 


H12 


3750.1 


+41 


3.03 + 0.13 












Hll 


3770.6 


+72 


5.01 + 0.16 












H10 


3797.9 


+108 


6.59 + 0.17 


6.55 










H9 


3835.4 


+56 


7.10 + 0.20 


7.16 










H8 


3889.1 


+59 


6.73 + 0.15 


7.00 










Call 


3933.7 


+191 


2.43 ± 0.09 


2.45 








IS 


He 


3970.1 


+70 


8.67 + 0.16 


8.56 










Hi 


4101.8 


+40 


8.17 + 0.24 


7.97 










H~) 


4340.5 




6.28 + 0.17 


6.02 












4O4U.0 














Nebular 


TT T i TT1 T TT 

Hel+leljll 


4385 




0.43 ± 0.07 


0.14 










He I 


4471.5 


+57 


0.40 ± 0.07 


1.09 










Mgll 


4481.1 


+77 


0.34 ± 0.08 










H/3 


4861.3 




6.55 ± 0.14 


6.08 










H/3 


4861.3 


+199 












Nebular 


He I 


4921.9 


+27 


0.56 ± 0.08 


0.35 










[OIH] 


4958.9 


+181 


-0.36 ± 0.06 










Nebular 


[OIII] 


5006.8 


+164 


-0.85 ± 0.05 










Nebular 


Mgl 


5170 




0.76 ± 0.12 


1.20 2 










Mgl 


5185.5 


+23 


0.43 ± 0.10 












Fe I+Ca I 


5270 




0.40 ± 0.04 


0.57 2 










Cr I+Cul 


5782 




0.75 + 0.08 


0.32 




1.09 ± 0.17 


2.53 




Nal 


5890.0 


+193 


6.20 + 0.10 


6.49 


+200 


6.59 ± 0.18 


6.74 


IS 


Nal 


5895.9 


+183 






+198 






IS 




6283.9 


+206 


1.52 + 0.11 




+174 


2.17 + 0.22 




DIB? 


[OI] 


6300.3 


+181 


-0.48 ± 0.06 










Nebular 




6497 




U. I £7 _1_ \J.\J t 


0.52 




n «o + 1 6 


0.71 




[Nil] 


6548.0 


+111 


-1.34 + 0.05 




+80 


-0.86 ± 0.08 




Nebular 


Ha 


6562.8 


+104 


-4.34 + 0.05 




+77 


-6.05 ± 0.08 




Nebular 


[Nil] 


6583.4 


+108 


-3.59 + 0.05 




+61 


-2.98 ± 0.09 




Nebular 


[SII] 


6716.5 


+126 


-1.34 + 0.06 




+65 


-1.57 + 0.09 




Nebular 


[SII] 


6730.9 


+121 


-1.10 + 0.06 




+61 


-1.14 + 0.09 




Nebular 


Til+Fel 


8468 




0.54 ± 0.08 






0.60 ± 0.08 






Call 


8498.0 


+31 


1.43 ± 0.13 3 




+133 


1.21 + 0.14 3 






Fel+Til 


8516 




0.41 ± 0.08 






0.83 ± 0.09 






Call 


8542.1 


+36 


4.09 ± 0.13 3 




+126 


4.14 + 0.14 3 






Pa 14 


8598.4 




1.51 + 0.16 












Call 


8662.1 


+38 


4.41 ± 0.13 3 




+135 


4.75 + 0.14 3 






Fel 


8688.6 


+38 


0.36 + 0.06 




+137 


0.51 ± 0.07 






Pa 12 


8750.5 




2.56 + 0.16 


2.45 




2.16 ± 0.18 


2.30 





1 all windows as denned by Bica & Alloin (1986) except for the following: 2 windows as defined by Faber et al. (1985); 
3 equivalcnt widths measured using the bands defined by Diaz et al. (1989). 



al. 1997). Typical Galactic paths with E(B-V)^1 produce 
central absorption depths of « 0.75 in the A6283 DIB (e.g. 
Jenniskens & Desert 1993, Herbig 1995); slightly stronger 
than what we observe in either cluster F (E(B— V)«l— 1.5) 
or cluster L (E(B-V)>2) in M82. The modest strengths of 
the DIBs along paths through the M82 ISM which contains 
strong Na I D lines and thus neutral gas suggests that we 
are observing these two clusters through regions with few 
ionization boundaries. In this case we might expect small 
scale structure to be common in the ISM, which could con- 
tribute to non-uniform interstellar obscuration across the 
faces of the clusters. Follow-up observations are needed to 
establish whether M82 is the second galaxy outside of the 



Local Group where DIBs have been detected (the first being 
supernova 1986G in NGC 5128; see Herbig 1995). 



4 COMPARISONS WITH THEORETICAL 
STAR CLUSTER MODELS 

4.1 Blue wavelength region 

Quantitative models now exist for spectra of single-age, solar 
metallicity star clusters which formed over short time scales. 
We used the results from the spectral synthesis code 'Projet 
d'Etude des Galaxies par Synthese Evolutive' (PEGASE; 
Fioc & Rocca-Volmerange 1997). These models are well- 
suited to our needs because they include supergiant stars 
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Figure 4. This plot shows the model spectrum for a star cluster with an age of 60 Myr computed with PEGASE. The flux is in units 



of erg 



per solar mass of cluster where a Salpeter initial mass function was assumed. The spectrum was synthesized using the 



Jacoby et al. (1984) library of stellar spectra and Geneva stellar evolution tracks 



that are critical in the spectra of young star clusters; mod- 
els that exclude such stars will not give proper matches to 
integrated spectra of younger super star clusters (e.g. Brodie 
et al. 1998). We selected the 'Geneva' stellar evolution tracks 
from Schaller et al. (1992) and Charbonnel et al. (1996) in 
computing cluster models with PEGASE. 

PEGASE also offers a choice of spectra synthesized from 
relatively low spectral resolution model atmospheres, or the 
Jacoby et al. (1984) spectrophotometry library. We chose 
the latter option for this project, as the near solar metallic- 
ity should be a reasonable match to young objects in M82, 
and we desired the best possible spectral resolution. Even 
though this appears to still be the highest spectral reso- 
lution Galactic spectrophotometric library that includes a 
wide range of spectral and luminosity classes, it has lower 
resolution (4.5A), and often lower per pixel signal-to-noise 
ratios, than our WHT spectra of M82 cluster F. Improved 
spectral libraries will be important for future comparisons 
between model and observed spectra of nearby SSCs. 

Before comparing with blue- region spectra of model star 
clusters, we smoothed the cluster F spectrum with a Gaus- 
sian having a = 1.0 pixel over AA3510-5002 A. This slightly 
degrades our resolution to about 2 A. We also shifted the 
observed spectrum by —175 km s _1 to account for the ap- 
proximate redshift of the stellar absorption lines. The IRAF 
package 'sarith' was then used to make ratios between the 
wavelength- and flux-calibrated data and lower spectral res- 
olution blue region PEGASE models. Models were computed 
for ages of 20, 40, 60, 80, and 100 Myr with a Salpeter IMF. 



The best-fitting 60 Myr model blue region spectrum is shown 
in Fig. | 

The comparisons between PEGASE models and data 
are very good. Ratios between observations and model spec- 
tra reveal both components of the interstellar H and K Ca II 
absorption lines, as well as fainter Balmer emission from the 
H7 and H<5 lines. Some 'ringing' is present among the higher 
Balmer series; probably due to the modest spectral resolu- 
tion of the spectral library as well as mismatches between 
wavelength scales. 

The ratios of spectra display an upward slope towards 
the red caused by dust obscuration in M82. OM78 estimated 
the reddening towards cluster F to be E(B— V)«l; our model 
fits suggest a slightly higher value of E(B— V)«1.5. However, 
a standard Galactic extinction model does not properly cor- 
rect the data; this law has too much curvature in the near 
UV; there is likely to be dust near cluster F and possibly 
even patchy obscuration (cf. conditions in the starburst nu- 
cleus of NGC 253; Watson et al. 1996), so it is not surprising 
that the Galactic extinction model fails (see also Calzetti 
1998). Because of this uncertainty in the proper choice of an 
obscuration correction, we cannot use the continuum spec- 
tral energy distribution at optical wavelengths to determine 
the cluster's age. 

Our efforts therefore concentrated on matching the 
Balmer absorption lines and gravity/temperature-sensitive 
Balmer jump as the basis for selecting best-fitting mod- 
els. While the models do not perfectly reproduce our ob- 
servations, they do not display gross difficulties of the type 
encountered by Brodie et al. (1998) in their efforts to fit 



© 1998 RAS, MNRAS 000, [I|-|h] 



8 J. S. Gallagher and L. J. Smith 




3600 3800 4000 4200 4400 

Wavelength (A) 



4600 



4800 



5000 



Figure 5. The ratio is plotted between the observed cluster F blue spectrum and 60 Myr PEGASE model displayed in Fig. |i| The ratio 
gives a model-dependent estimate for the mass of the cluster required to provide the observed flux, and in the absence of dust obscuration 
should be a constant function of wavelength. This ratio is normalized such that a ratio of 1 would correspond to a model cluster mass 
of 1.6xl0 6 M Q 



blue spectra of SSCs in NGC 1275. We therefore had no 
need to invoke an unusual form of the initial mass func- 
tion; this cluster is fitted adequately with a Salpeter IMF. 
The 40 Myr PEGASE model gives the best fit to the H/3, 
H7, and H<5 absorption lines but does not do as well for the 
Balmer jump region as the 60 and 80 Myr models. The 20 
and 100 Myr cluster models gave worse fits in all areas. We 
therefore adopt an age for cluster F of 60±20 Myr from the 
PEGASE model fitting. 

The ratio between the observed cluster F spectrum and 
60 Myr model is shown in Fig. |^. Since we are comparing a 
monochromatic luminosity with a predicted luminosity per 
unit cluster mass, the ratio is in units of solar masses. The 
normalization leads to a ratio of 1 in the plot correspond- 
ing to a cluster mass of 1.6 x 10 6 Mq . If we correct by a 
factor of 40-200 for As =4-6 magnitudes, then we see that 
a model cluster mass of >10 7 Mq is needed to reproduce 
the blue luminosity. Because the form of the initial mass 
function in SSCs is poorly known, this mass is highly uncer- 
tain. However, it does indicate that cluster F is likely to be 
very massive, in agreement with results from earlier studies 
(OM78, O'Connell et al. 1995). 



4.2 Red wavelength region 

The prime diagnostic for measuring ages in the red part of 
the spectrum is the strength of the Call triplet (CaT). The 
presence of a strong CaT is often used as an indicator of 
red supergiants (RSG) and therefore an age of 10 7 yr. The 



CaT is, however, strong in yellow supergiants as well, and as 
discussed by Mayya (1997), is not an unambiguous signature 
of RSGs. Mayya (1997) finds that, at solar metallicity, the 
evolution of the strength of the two strongest members of 
the CaT (Ca2+Ca3) shows two peaks - one from RSGs at 
~ 10 7 yr (W\ ~ 11 A), and one at 60 Myr due to red giants 
and asymptotic branch stars (W\ w 7 A). 

In M82 clusters F and L, the combined equivalent 
widths of Ca2+Ca3 are 8.5 and 8.9 A respectively, uncor- 
rected for Paschen line absorption. Spectral synthesis mod- 
els (Leitherer, priv. comm.) at ages of 40-80 Myr predict 
equivalent widths of 7.9-8.3 A (including the Paschen lines) 
and 5.6-5.9 A (pure CaT). The strength of the CaT in the 
spectrum of SSC F is therefore in accord with the age de- 
duced from the blue spectrum of ~ 60 Myr. For SSC L, the 
only potential age information that we have is the strength 
of the CaT. The similarity in the strength of this feature, 
and the overall spectral appearance, indicates that SSC L is 
probably of a similar age to F. 



5 DISCUSSION 

The combination of high quality optical spectra and good 
spectral synthesis models for star clusters has allowed us to 
derive an age for M82 SSC F. Since SSCs are hallmarks of 
episodes of intense star formation (see Hunter, O'Connell, & 
Gallagher 1994; O'Connell et al. 1994; Watson et al. 1996), 
cluster F provides one benchmark in the history of the M82 
starburst. While this is only one star cluster, it is a massive 
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object, and so its birth would likely have been noticeable on 
a galactic scale; the SFR for making this cluster during a 
1 Myr dynamical time scale would have been >1 Mq yr _1 . 
Thus our data suggest a minimum age for the starburst in 
M82 of 40 Myr and a most probable age of > 60 Myr. 

Ages of the M82 starburst have been previously derived 
from integrated spectra of the central regions of this galaxy 
by Rieke et al. (1993). Their preferred models consist of two 
bursts within the past ~30 Myr, each with a duration of 
about 5 Myr. Doane & Mathews (1993) based their models 
on the energetics of the stellar populations, including su- 
pernova rates estimated from radio observations, and found 
the starburst has been active for 30-60 Myr. Satyapal et 
al. (1997) observed twelve luminous concentrations within 
the central starburst in the infrared. They identify these ob- 
jects with embedded SSCs, whose frequent presence is also 
required to fit the numbers of radio supernovae remnants 
(Golla, Allen, & Kronberg 1996). The ages of these SSCs 
with high near-IR luminosities are estimated by Satyapal 
et al. (1997) to be near 10 Myr, and they present a model 
in which the starburst has propagated outwards from the 
nucleus during this time period. 

At larger radii there is evidence that the starburst is 
older. For example, Shen & Lo (1995) measured a change in 
the properties of the molecular medium between the central 
and mid-disk zones of M82. The molecular gas at a radius 
of 390 pc is much more disturbed than that in the inner 
125 pc region, suggesting that the outer gas was disrupted 
by an earlier starburst. They identify this with the 30 Myr 
starburst event suggested by Rieke et al. (1993). OM78 also 
noted that the outer region is probably an older event than 
that in the centre. 

The structure of M82 is an important factor when dis- 
cussing the "age" of the starburst (see Telesco 1988). Cluster 
F lies to the southwest of the M82 nucleus, slightly above its 
very heavily obscured mid-plane (Larkin et al. 1994). The 
rapid variation in extinction across this area is also illus- 
trated by the presence of the highly-reddened SSC L near 
cluster F. It is in the region of M82 where the near-IR surface 
brightness remains roughly constant, and may be associated 
with a ring structure that is inferred to surround a central 
bar (see Telesco et al. 1991, Larkin et al. 1994, Achtermann 
& Lacey 1995). 

We can also use the measured radial velocities of the 
absorption and emission lines in the spectra of clusters F 
and L to obtain information on their locations within M82. 
Gotz et al. (1990) and McKeith et al. (1993) present detailed 
studies of the structure of the disk of M82 using optical and 
near-IR emission and absorption line velocities. In common 
with previous studies, they find that the size of the rotational 
velocity gradient increases with wavelength because of the 
obscuration in M82. The heliocentric velocities we derive for 
the interstellar absorption lines (+189 ± 5kms _1 ) and the 
blue nebular emission lines (+181 ± 12kms _1 ; Sect. 3.3) 
agree well with those determined for the same lines near the 
position of cluster F by Gotz et al. (1990). These velocities 
are below the systemic value, but well above the rotation 
velocity in this region of M82 (e.g. Weliachew et al. 1984, 
Achterman & Lacey 1995), consistent with a model where 
the emission lines originate in ionized gas located on the 
near side of the outer disk of M82 (McKeith et al. 1993). 

In contrast, the radial velocity of +35 km s _1 we derive 



for cluster F is very near the expected rotation curve peak 
at its position. Cluster F therefore should be located near 
the mid-point of our sight line through the disk of M82, at 
a radial distance of about 440 pc from the nucleus, near the 
edge of the dense gaseous ring. Possibly a hole in the dust 
layer is allowing us to see cluster F relatively deep within 
M82 (see Larkin et al. 1994). For cluster L, we find a CaT 
velocity of +131kms _1 which is similar to our red HII re- 
gion emission line velocities and those of Gotz et al. (1990). 
This suggests that L is located towards the outer edges of 
the disk at a radius of somewhat more than 440 pc in a 
region of high obscuration. 

The different ages derived for the M82 starburst then 
may simply reflect variations of the star formation history 
within this extraordinary galaxy. If we assume that cluster 
F is a typical component of the 'ring' or mid-disk, then the 
starburst occurred in this region rather long ago, and in the 
last 20-30 Myr has become more concentrated in the central 
part of the galaxy, where a strong bar may be present. This 
model is consistent with previous results in suggesting the 
earlier M82 starburst occurred in the mid-disk, and the more 
recent event is more centrally concentrated. Age dating of 
a larger population of star clusters beyond the currently 
active central-region starburst has the potential to provide 
a clearer view of the development of the M82 starburst. 



6 CONCLUSIONS 

We have obtained high signal-to-noise ratio optical spectra 
with good wavelength resolution for star clusters F and L 
in M82. Quantitative fits to theoretical models are success- 
ful in reproducing the spectra with an SSP having an age of 
~ 60 Myr for cluster F. This corresponds to a main sequence 
turnoff mass of about 7 Mq, consistent with the tempera- 
tures of stars associated with weaker spectral lines seen in 
the blue region. Cluster L appears to have a similar age, 
but the result is much less certain due to the lack of blue 
spectra. The combination of high quality spectra and stellar 
population models which include supergiant stars provides 
a powerful means to analyze young super star clusters. 

Combining our estimate of V=15.8 with an extinction 
of E(B— V)^ 1.5 we derive My = —16.5 for cluster F, two 
magnitudes brighter than that determined by O'Connell et 
al. (1995). It therefore appears likely that these clusters are 
extremely massive and may have the ability to survive as 
bound systems over long time scales (e.g. Goodwin 1997). A 
measurement of the stellar velocity dispersion in cluster F 
would be extremely useful in constraining its mass and thus 
similarity to globular star clusters as well as future evolution 
(cf. Ho & Filippenko 1996). 

The spectra indicate that while cluster L is more ob- 
scured than cluster F, the latter star cluster likely lies near 
the outer edge of the main starburst region in M82. It 
is an extremely luminous super star cluster, whose pres- 
ence strongly suggests that the M82 starburst was quite ac- 
tive 60 Myr ago. These clusters also offer useful test paths 
through the ISM of M82, and in addition to the well-known 
atomic absorption of Nal D and Call H+K, we find a prob- 
able DIB feature. 
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